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Sphere Packings II 


Thomas C. Hales 


Abstract: An earlier paper describes a program to prove the Kepler conjecture on 
sphere packings. This paper carries out the second step of that program. A sphere 
packing leads to a decomposition of into polyhedra. The polyhedra are divided 
into two classes. The first class of polyhedra, called quasi-regular tetrahedra, have 
density at most that of a regular tetrahedron. The polyhedra in the remaining class 
have density at most that of a regular octahedron (about 0.7209). 

Section 1. Introduction 

This paper is a continuation of the first part of this series [4]. The terminology and 
notation of this paper are consistent with this earlier paper, and we refer to results 
from that paper by prefixing the relevant section numbers with ‘T. 

We review some definitions from [4]. Begin with a packing of nonoverlapping spheres 
of radius 1 in Euclidean three-space. The density of a packing is defined in [1]. It 
is defined as a limit of the ratio of the volume of the unit balls in a large region 
of space to the volume of the large region. The density of the packing may be 
improved by adding spheres until there is no further room to do so. The resulting 
packing is said to be saturated. 

Every saturated packing gives rise to a decomposition of space into simplices called 
the Delaunay decomposition [8]. The vertices of each Delaunay simplex are centers 
of spheres of the packing. By the definition of the decomposition, none of the centers 
of the spheres of the packing lie in the interior of the circumscribing sphere of any 
Delaunay simplex. We refer to the centers of the packing as vertices. Vertices that 
come within 2.51 of each other are called close neighbors. 

The Delaunay decomposition is dual to the well-known Voronoi decomposition. If 
the vertices of the Delaunay simplices are in nondegenerate position, two vertices 
are joined by an edge exactly when the two corresponding Voronoi cells share a 
face, three vertices form a face exactly when the three Voronoi cells share an edge, 
and four vertices form a simplex exactly when the four corresponding Voronoi cells 
share a vertex. In other words, two vertices are joined by an edge if they lie on a 
sphere that does not contain any other of the vertices, and so forth (again assuming 
the vertices to be in nondegenerate position). 

We say that the convex hull of four vertices is a quasi-regular tetrahedron (or simply 
a tetrahedron) if all four vertices are close neighbors of one another. If the largest 
circumradius of the faces of a Delaunay simplex is at most v^, we say that the 
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simplex is small Suppose that we have a conhguration of six vertices in bijection 
with the vertices of an octahedron with the property that two vertices are close 
neighbors if and only if the corresponding vertices of the octahedron are adjacent. 
Suppose further that there is a unique diagonal of length at most 2^2. In this 
case we call the convex hull of the six vertices a quasi-regular octahedron (or simply 
an octahedron). A Delaunay star is dehned as the collection of all quasi-regular 
tetrahedra, octahedra, and Delaunay simplices that share a common vertex v. 

We assume that every simplex S in this paper comes with a hxed order on its 
edges, 1,... ,6. The order on the edges is to be arranged so that the hrst, second, 
and third edges meet at a vertex. We may also assume that the edges numbered 
i and i 3 are opposite edges for i = 1, 2, 3. We dehne S{yi,... , y^) to be the 
(ordered) simplex whose ith edge has length If S' is a Delaunay simplex in a 
hxed Delaunay star, then it has a distinguished vertex, the vertex common to all 
simplices in the star. In this situation, we assume that the edges are numbered so 
that the hrst, second, and third edges meet at the distinguished vertex. 

A function, known as the compression r(S), is dehned on the space of all Delaunay 
simplices. Set hoct = (—Stt-I- 12 arccos(l/v^))/v^ ~ 0.720903. Let S'be a Delaunay 
simplex. Let B be the union of four unit balls placed at each of the vertices of S. 
Dehne the compression as 

r(S) = —5octvol(S) -|- vol(S n B). 

We extend the dehnition of compression to Delaunay stars D* by setting r(D*) = 
^ r(S), with the sum running over all the Delaunay simplices in the star. We dehne 
a point (abbreviated pt) to be r(S(2, 2, 2, 2, 2, 2)) 0.0553736. The compression is 

often expressed as a multiple of pt. 

There are several other functions of a Delaunay simplex that will be used. The 
dihedral angle dih(S) is dehned to be the dihedral angle of the simplex S along 
the hrst edge (with respect to the hxed order on the edges of S). The solid angle 
(measured in steradians) at the vertex joining the hrst, second, and third edges is 
denoted so^A). Let rad(5') be the circumradius of the simplex S. More generally, 
let rad(F) denote the circumradius of the face of a simplex. Let ? 7 (a, b, c) denote 
the circumradius of a triangle with edges a, 6, c. Explicit formulas for all these 
functions appear in 1.8. 


Fix a Delaunay star D* about a vertex vq, which we take to be the origin, and we 
consider the unit sphere at vq. Let vi and V 2 be vertices of D* such that uq, fi, 
and V 2 are all close neighbors of one another. We take the radial projections pi of 
Vi to the unit sphere with center at the origin and connect the points pi and p 2 
by a geodesic arc on the sphere. We mark all such arcs on the unit sphere. The 
closures of the connected components of the complement of these arcs are regions 
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on the unit sphere, called the standard regions. We may remove the arcs that do 
not bound one of the regions. The resulting system of edges and regions will be 
referred to as the standard decomposition of the unit sphere. 

Let C be the cone with vertex vq over one of the standard regions. The collection 
of the Delaunay simplices, quasi-regular tetrahedra, and quasi-regular octahedra 
of D* in C (together with the distinguished vertex uq) will be called a standard 
cluster. Each Delaunay simplex in D* belongs to a unique standard cluster. 

A real number, called the score, will be attached to each cluster. Each star receives 
a score by summing the scores for the clusters in the star. 

The steps of the Kepler conjecture, as outlined in Part I, are 

1. A proof that even if all standard regions are triangular, the total score is less 
than 8pt 

2. A proof that the standard clusters with more than three sides score at most Opt 

3. A proof that if all of the standard regions are triangles or quadrilaterals, then 

the total score is less than 8 pt (excluding the case of pentagonal prisms) 

4. A proof that if some standard region has more than four sides, then the star 

scores less than 8 pt 

5. A proof that pentagonal prisms score less than 8pt 

The proof of the hrst step is complete. The other steps are briefly discussed in Part 
I. This paper establishes step 2. Partial results have been obtained for step 3 [5]. 
C.A. Rogers has shown that the density of a regular tetrahedron is a bound on the 
density of packings in [8]. The main result of this paper may be interpreted as 
saying that the density {doct ~ 0.7209) of a regular octahedron is a bound on the 
density of the complement in of the quasi-regular tetrahedra in the packing. 

The score of a Delaunay star is obtained by mixing Delaunay stars with the dual 
Voronoi cells. Delaunay stars D* and the associated function P behave much better 
than estimates of density by Voronoi cells, provided each Delaunay simplex in the 
Delaunay star has a small circumradius. Unfortunately, r(5') gives an increasingly 
poor bound on the density as the circumradius of the Delaunay simplex S increases. 
When the circumradius of S is greater than about 1.8, it becomes extremely difficult 
to prove anything about sphere packings with the function r(>S'). The score is 
introduced to regularize the irregular behavior of r(5'). 

Voronoi cells also present enormous difficulties. The dodecahedron shows that a 
single Voronoi cell cannot lead to a bound on the density of packings better than 
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about 0.75. This led L. Fejes Toth to propose an approach to the Kepler conjecture 
in which two layers of Voronoi cells are considered: one central Voronoi cell and 
a number of surrounding ones. Wu-Yi Hsiang has made some progress in this 
direction, but there remain many technical difficulties [3], [7]. 

The method of scoring in this paper seeks to combine the best aspects of both 
approaches. When the circumradius of a simplex is small, we proceed as in Part I. 
However, when the circumradius of a simplex is large, we switch to Voronoi cells. 
Remarkably, these two approaches may be coherently combined to give a meaningful 
score to Delaunay stars and, by extension, a bound on the density of a packing. 
The calculations of this paper suggest that this hybrid approach to packings retains 
the best features of both methods with no (foreseeable) negative consequences. 

Section 2. Some polyhedra 


Sometimes the tip of a Voronoi cell protrudes beyond the face of a corresponding 
Delaunay simplex (see Diagram 2.1.a). This section describes a construction that 
amounts to slicing off the protruding tip of a Voronoi cell and reapportioning it 
among the neighboring cells (see Diagram 2.1.b). 



Let D* be a Delaunay star with center vq = 0. Let V be the Voronoi cell around 
uq, obtained by duality from D*. As a matter of convenience, we may assume that 
each point in belongs to a unique Voronoi cell by making an arbitrary choice for 
each point on the boundary of a cell. If R is a standard cluster (possibly a single 
quasi-regular tetrahedron) in D*, let C{R) denote its cone over vq: 

C{R) = {tx : t > 0, X E R}. 


In general, V fl C{R) depends on more vertices than just those in the cluster R. 
It is convenient to consider the slightly larger polyhedron dehned by just the 
vertices of D* that are in R. That is, let be the intersection of C(R) with the 
half-spaces {x : x ■ Vi < Vi ■ Vijl^ 'ii ^ 0}, where {vi}i are the vertices (other than 
un) of the simplices and quasi-regular solids in the cluster R. The faces of VS at vn 
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are contained in the triangnlar faces bonnding the standard region of R. The other 
faces of are contained in planes throngh the faces of the Voronoi cell V. We refer 
to these as Voronoi faces. If R is not a qnasi-regnlar tetrahedron, set Vr = If i? 
is a qnasi-regnlar tetrahedron, we take the slightly smaller polyhedron Vr obtained 
by intersecting with the half-space (containing no) bonnded by the hyperplane 
throngh the face of R opposite the origin vq- (This may cnt a tip from the Voronoi 
cell.) By constrnction, Vr depends only on the simplices in R. The polyhedron Vr 
is based at the center of some Delannay star, giving it a distingnished vertex v. We 
write Vr = VR(n) when we wish to make this dependence explicit. 

By constrnction, V^ D V fl C(R). It is often trne that Vr = V fl C(R). Let ns 
stndy the conditions nnder which this can fail. We say that a vertex w clips a 
standard clnster R (based at ni) if w vi and some point of VR(ni) belongs to the 
Voronoi cell at w. Part I makes a thorongh investigation of the geometry when a 
vertex w clips a qnasi-regnlar tetrahedron. (The vertex w mnst belong to a second 
qnasi-regnlar tetrahedron that shares a face with S. The shared face mnst have 
circnmradins greater than and so forth.) 

Lemma 2.2. Let R, based at a vertex no, be a standard clnster other than a 
qnasi-regnlar tetrahedron. Snppose it is clipped by a vertex w. Then there is a face 
(no, vi,V 2 ) of R snch that (rc, no, ni, n 2 ) is a qnasi-regnlar tetrahedron. Fnrthermore, 
(no,ni,n 2 ) is the nniqne face of the qnasi-regnlar tetrahedron of circnmradins at 
least \/2. 


Proof: Consider a point p in Vr \ V. Then there exists a vertex w ^ C{R) of D* 
snch that p ■ w > w ■ w/2. The line segment from p to m intersects the cone C{F) 
of some triangnlar face F that bonnds the standard region of R and has no as a 
vertex. Let ni and V 2 be the other vertices of F. By the constrnction of the faces 
bonnding a standard region, the edges of F have lengths between 2 and 2.51. 

Consider the region X containing p and bonnded by the planes Hi = span(ni,rc), 
H 2 = span(n 2 , m), iLs = span(ni, n 2 ), = {a; : a; ■ ni = ni ■ ni/2}, and = {a: : 

X ■ V 2 — V2 ■ V2/ 2 }. The planes 1^4 and contain the faces of the Voronoi cell at 
no dehned by the vertices vi and V 2 - The plane iLs contains the face F. The planes 
Hi and H 2 bonnd the region containing points, snch as p, that can be connected 
to m by a segment that passes throngh C{F). 

Let P = {x : X ■ w > w ■ w/2}. The choice of w implies that V fl P is nonempty. 
We leave it as an exercise to check that V fl P is bonnded. If the intersection 
of a bonnded polyhedron with a half-space is nonempty, then some vertex of the 
polyhedron lies in the half-space. So some vertex of X lies in P. 

We claim that the vertex of X lying in P cannot he on Hi. To see this, pick 
coordinates (xi , xo) on the plane Pi with origin nn = 0 so that ni = (0, z) (with 
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z > 0) and X D Hi C X' := {{xi^X 2 ) : iCi > 0, X 2 < z/2}. See Diagram 2.3. If 
X' meets P, then the point vi/2 lies in P. This is impossible, becanse every point 
between no and vi lies in the Voronoi cell at no or ni, and not in the Voronoi cell 
of w. (Recall that \vi — no| < 2.51 < 2v^.) 



Similarly, the vertex of X in P cannot lie on H 2 - Thns, the vertex mnst be the 
nniqne vertex of X that is not on Hi or H 2 , namely, the point of intersection of 
P 3 , H 4 , and H 5 . This point is the circnmcenter c of the face F. We conclnde 
that the polyhedron Xq := X D P contains c. Since c G Xq, the hypotheses of 
Lemma 1.3.4 are met for T = F, and the vertices no, ni, V 2 , and w are the vertices 
of a qnasi-regnlar tetrahedron S. By 1.3.4, \w' — Vi\ < 2.3, for i = 0,1,2. The 
circnmradins of the face F is between and 2.51/v^ ~ 1.449. I 

In the same context, if w and w' both clip R, then the regions they cnt from 
Vnivo) are disjoint. For otherwise, a common point wonld belong to both Vg{w) 
and Vg{w') where S and S' are the two qnasi-regnlar tetrahedra constrncted by 
the Lemma. Part 1.3 shows that S and S' share their nniqne face of circnmradins 
greater than \/ 2 . This is impossible, becanse the le mm a states that this face is 
shared with R. 

Althongh the polyhedron Xq belongs to the Voronoi cell at w, it is inclnded in the 
polyhedron Vr. Similarly, by repeating the constrnction at vi and ^ 2 , we hnd that 
there are small regions Vi, X 2 (with vertex c) in polyhedra and Vr^ at vi and 
V 2 respectively that belong to the Voronoi cell at w. 

Call the nnion V 0 UV 1 UV 2 the tip protrnding from the qnasi-regnlar tetrahedron S. 
Associated with a qnasi-regnlar tetrahedron is at most one such tip. (The tip must 
protrude from the face of S with circnmradins greater than ^/2.) By construction, 
the tip is the set of points 


{x : \x — w\ <\x — Vi\^ fori = 0 , 1 , 2 ; det(a:, ui, ^ 2 ) det(r(;, ui, ^ 2 ) < 0 }. 


This is Vl{w)\Vs{w). 
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The tip is a subset of the Voronoi cell at w. Section 1.3 explains the conditions 
under which this can fail to hold. There must be another vertex u ^ w with the 
property that \u — Vi\ < 2.3, for f = 0,1, 2. Then u, vq, vi, and V 2 are the vertices 
of a second quasi-regular tetrahedron S' with face F, and this is contrary to our 
assumption that R is not a quasi-regular tetrahedron. 


Corollary 2.4. The polyhedra cover evenly as we range over all the standard 
clusters of all the Delaunay stars of the packing. 


Proof: The preceding analysis shows that the polyhedra Vr are obtained from the 
Voronoi cells by taking each protruding tip, breaking it into three pieces Xq, Vi, 
X 2 , and attaching the piece Xi to the Voronoi cell at Vi. The Voronoi cells cover 
evenly. As a result of this analysis, we see that the polyhedra Vji cover evenly.B 


To give one example of the size of the tip, we consider the extreme case of the 
tetrahedron S = 5'(2, 2, 2, 2.51, 2.51, 2.51). Diagram 2.5 shows a correctly scaled 
drawing of a tip protruding from the largest face of S. 



Section 3. The score attached to a Delaunay star 


This section gives some rules for computing the score. They were developed as 
a result of computer experimentation suggesting when it is advantageous to use 
Voronoi cells over Delaunay simplices. This section actually gives an entire family 
of scoring systems. This extra bit of flexibility will be useful as we encounter new 
examples in the remaining steps of the program. We expect the score to satisfy 
Conjecture 1.2.2, which asserts that the score of a Delaunay star is at most Spt, 
for all the scoring systems satisfying properties 1-4 below. The Kepler conjecture 
is true if Coniecture 1.2.2 holds for anv one such scoring svstem. We have found 





through experimentation that small or seemingly innocent changes in the score can 
lead to enormous changes in the complexity of the optimization problem. 

3.1. This paper proves the second step of the program for all of the scoring systems 
presented below. Write (j{S) for the score of S. 

1. Suppose that the standard cluster i? is a single quasi-regular tetrahedron: R = S. 
When the circumcenter of S is contained in S, 

-4hoctvol(Vi? n C{S)) -s 4sol(S')/3 

is an analytic f un ction of the lengths of the edges. This expression has an analytic 
continuation, denoted vor(S', Vr), to simplices S that do not necessarily contain 
their circumcenter. 

If rad(S') > 1.41, then dehne the score to be vor(S', Vr). If rad(S') < 1.41, then 
dehne the score to be the compression r(S'). (This rule agrees with the dehnition 
of vor(S') given in Section 1.2.) 

2. Let S' be a small simplex that is not a quasi-regular tetrahedron. The score of 
S will be either vor(S) or r(S) depending on criteria to be determined by future 
research.* These criteria may depend on whether S belongs to a quasi-regular oc¬ 
tahedron, but not on the position of any vertices of the packing outside S. It is 
essential for the scoring at all four vertices to have the same type (Voronoi or com¬ 
pression) . The only constraint imposed by the second step of the Kepler conjecture 
will be cr(S) < 0, if S is small. This leads to the following mild restrictions on the 
use of Voronoi scoring. 

If one of the hrst three edges is the long edge (say the hrst), compression scoring is 
to be used if the second, third, and fourth edges have length at most 2.06, and the 
hfth and sixth edges have length at most 2.08. 

If one of the last three edges (say the fourth) is the long edge, compression scoring 
is to be used if (a), (b), (c), and (d) hold. 

(a) The hrst edge has length at most at most 2.06. 

(b) The second and third edges have length at most 2.08. 

(c) The hfth and sixth edges have length at most 2.2. 

(d) The fourth edge has length at most 2.58, or the hfth and sixth edges have 
lengths at most 2.12. 

*More generally, we might add a small constant c to the score of S at one of its vertices and 
subtract the same constant from another vertex. 



9 


3. Suppose that R is any standard cluster other than a quasi-regular tetrahedron. 
The cluster is a union of Delaunay simplices Si^... , Sr- Index the simplices so that 
S'!,... , S'p, for some p < r are the small simplices in the cluster. We dehne the 
score of the cluster R to be 

a{Si)+ Y voY{Si,VR), 

l<z<p p<.iKr 

where vor(S', Vr) = 4(—hoctvol(VR fl C{S)) -|- sol(S')/3). 

4. If D* is a Delaunay star, then its total score u{D*) is a sum of the scores of the 
standard clusters of D*. 


Consider the quasi-regular tetrahedron S of Section 2 with vertices uq, ui, ^ 2 , and 
w that has a protruding tip Xq U U X 2 - Let sol^;(S') denote the solid angle of S 
at the vertex v. The analytic continuation vor(S', Vr) has the following geometric 
interpretation. 

vor(S', Vr{v)) = -45oct(vol(S', -F A{v)) + 4 soC(S')/3, 

with the correction term A{vi) = —vol(W )5 for * = 0,1, 2, and 

3 

A{w) = y^vol(W)- 

i=l 


The only pieces that are compression scored are small simplices, everything else 
is Voronoi scored. The small simplices that are compression scored will be called 
simplices of compression type. The Voronoi-scored small simplices will be called 
simplices of Voronoi type. We dehne the restricted cell of a cluster R to be the 
complement in Vr of the small simplices in the packing. 


Lemma 3.2. (1) The score of a cluster depends only on the cluster, and not on 
the way it sits in a Delaunay star or in the Delaunay decomposition of space. 

(2) Let A denote the vertices of a saturated packing. Let A^v denote the vertices 
inside the ball of radius N. (Fix any center for the ball.) Let D*{v) denote the 
Delaunay star at n G A. Then the score satishes (in Landau’s notation) 

5^ <7 (D’(v)) = 5^ r(D'(v)) + O(N^). 

An An 


T>_ r. _ 4. U,, __^4- 
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(2) The score reapportions the compression of a given star among snrronnding 
stars. The second part of the lemma follows from the claim that everything is 
acconnted for, if we ignore the bonndary effects cansed by the trnncation N. Space 
is partitioned into regions each connted —ASoct times by the compression of some 
star. Each point in a sphere of the packing is connted fonr times by the compression 
of some star. To verify 3.2.2, we mnst check that the same holds of the score. 

We switch from Voronoi to compression scoring on certain small simplices. The 
faces F of a small simplex S satisfy rad(F) < so no point on a face F of S 
can be closer to another vertex in the packing than it is to the closest vertex of F. 
This has two implications. First, the only polyhedra Vr meeting a small simplex S 
are the fonr based at the vertices of S. Second, Let F be a standard clnster. Let S 
be a small simplex in R. Then Vr fl C{S) = VrC] S. (In other words, tips cannot 
protrnde from a small simplex.) This means that the restricted cells and small 
simplices cover space evenly. This decomposition is compatible with the standard 
decomposition of a Delannay star. 

Consider the rnles dehning the score. In connting the part of the volnme of a sphere 
contained in a simplex S, we see that it appears fonr times with weight 1 for a total 
weight of 4, when S' is a small simplex of compression type. It appears once with 
weight 4 for a total weight of 4, when S is of Voronoi type. 

The resnlt is now clear. I 


Remark 3.3. It is nsefnl to su mm arize the proof from a slightly different point 
of view. If S' is a qnasi-regnlar tetrahedron or a small Delannay simplex, then 
the snm of its fonr scores, for each of its fonr vertices, is 4r(S). This follows 
directly from the dehnitions (and the proof of Lemma 3.2) if the circnmcenter of 
S is contained in S (which is always the case for small simplices), and it follows 
by analytic continnation in general. Any other point in space belongs to a nniqne 
Voronoi cell centered at some vertex v. If the point is not in a tip protrnding from 
a qnasi-regnlar tetrahedron, it is connted in the score at v. If, however, the point 
belongs to a protrnding tip, it is connted in the score at exactly one of the three 
vertices, other than v, of the qnasi-regnlar tetrahedron. In this way, every point in 
is acconnted for. 

Remark 3.4. The choice of the parameter n = 1.41 in Rnle 1 is somewhat arbi¬ 
trary. The choice is based on the comparison of the fnnctions 

f^(x) = vor(A(2, 2, 2, 2.51,2.51,0:), Vs) and f 2 ix) = r(A(2, 2, 2, 2.51, 2.51, x)). 

The difference fi{x) — f 2 {x) has a zero for some x G [2.2603,2.2604]. This gives a 
crnde estimate of when it is advantageons to switch from r(5') to vor(5'. Vs). The 
constant 1.41 is a little more than rad(5'(2, 2, 2, 2.51, 2.51, 2.2604)) ^ 1.405656. 



11 


Proposition 3.5. The Delaunay stars in the face-centered cubic and hexagonal- 
close packings score 8 pt. 

Proof: The eight regular tetrahedra each score 1 pt, and each regular octahedron 
scores Opt, because it has density Soct, for a total of 8 pt. I 

We will see in Proposition 4.6 that the regular octahedron can be broken into 
smaller pieces that score 0 pt. 

Section 4. The Main Theorem 


Theorem 4.1 (a) The score of any small quasi-regular tetrahedron is at most 1 pt. 
(b) The score of any other standard cluster is at most 0 pt. 

Proof. Statement (a) is a special case of Calculation 1.9.1. A quasi-regular tetra¬ 
hedron of Voronoi type scores less than 0 pt by Lemma 1.9.17. In the remainder of 
the proof, we actually prove a much stronger statement. We explicitly decompose 
each cluster (other than a quasi-regular tetrahedron) into a number of pieces and 
show that the density of each piece is at most doct- Since vor(S', Vn) and r(S') are 
zero precisely when the corresponding densities are 5oct (or when the volumes are 
zero), the theorem will follow. The relevant pieces will be congruent to one of the 
following types: 

1 . A small simplex that is not a quasi-regular tetrahedron 

2 . A set {tx : 0 < t < 1, a:GP 2 }cM^, where P 2 is a measurable set and every 
point of P2 has distance at least 1.18 from the origin (Diagram 4.2.a) 

3. A set {tx : 0 < t < 1, x E P3} C where P 3 is a wedge of a disk of the form 

P 3 = {{xi,X 2 , X 3 ) : X 3 = zq, x\ + xl< 2 , 0 <X 2 < axi), 
for some a > 0 and some 1 < zq < 1.18 (Diagram 4.2.b) 

4. A Rogers simplex R{a, b, a/2) where 1 < a < 1.18 and 4/3 < 6 ^ < 2 (see Section 
1.8.6 and Diagram 4.2.c) 
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Diagram 1^.2 





In the first type, a unit ball is placed at each vertex of the simplex S', and the 
density is the ratio of the volume of the part of the balls in S to the volume of S. 
In the second, third, and fourth types, a unit ball is placed at the origin, and the 
density is the ratio of the volume of the part of the ball in the region to the volume 
of the region. 

We decompose all of into these four types and quasi-regular tetrahedra. Set all 
the quasi-regular tetrahedra aside. Classify all the small simplices, including those 
contained in a quasi-regular octahedron, as regions of the hrst type. There remain 
the restricted cells. Now £x a Delaunay star D*, with center at the origin, and 
consider the restricted cell of one of the clusters in the star. We may assume that 
the restricted cell does not he in a quasi-regular tetrahedron. Break the restricted 
cell up further by taking its intersection with the cones over each of its Voronoi 
faces F. Let X be one such intersection. If the face F has distance more than 1.18 
from the origin, classify W as a region of the second type. Now assume the face 
F has distance h at most 1.18 from the center. Because h < v^, the point in the 
plane of F closest to the origin lies on the face F. The set of points P 2 on the face 
F at distance greater than from the origin gives rise to a region of the second 
type. To study what remains, we may truncate F by intersecting it with a ball of 
radius \/2. Let F' C F be the truncated face. 

By Voronoi-Delaunay duality, the face F' lies in the bisecting plane between 0 and 
some vertex v of the Delaunay star. Consider the collection of triangles formed by 
0, V, and another vertex of the Delaunay star D*, with the property that either 
the triangle has circumradius at most \/2 or all three edges of the triangle have 
lengths between 2 and 2.51. Consider the half-planes (bounded by the line through 
0 and v) containing the various triangles in this collection. This fan of half-planes 
partitions the face F' into a collection of wedge-shaped pieces. Consider one of 
them, F". We claim that it has the form of Diagram 4.3. 
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More precisely, F" is bounded by two triangular faces Fi and F 2 (in this collection 
of triangles), two edges ci and 62 of the Voronoi cell dual to the triangles, and an 
arc a obtained from the truncation. The two edges ei and 62 are perpendicular 
to the faces Fi and F 2 , respectively, by the dehnition of Voronoi-Delaunay duality. 
The edges e\ and 62 meet the faces Fi and F 2 , respectively, by the construction 
of restricted cells. The edges ei and 62 cannot intersect at any point less than 
from the origin, because the point of intersection would be a point equidistant from 
the four vertices of a simplex formed by the vertices of Fi and F 2 - The simplex 
would have circumradius less than \/2. Its faces would then also have circumradius 
less than \/2, so that the Delaunay simplex is small. This is impossible, since all 
small simplices have already been classihed as regions of the hrst type. 

Lemma 4.4. In this context, assume that the faces Fi and F 2 form an acute angle, 
and let p be the point at which the line through ei meets the plane through F 2 . Let 
wi and W 2 be the third vertices of the faces Fi and F 2 respectively (that is, those 
other than 0 and v). If the distance from 0 to p is at most v^, then the simplex 
( 0 , n, rci, ^ 2 ) is small or a quasi-regular tetrahedron. 

Proof: Suppose that the distance from p to 0 is at most v^. Let ci be the 
circumcenter of Fi. Let S be the simplex (0, n, wi, W 2 )- 

We claim that p lies in the interior of the triangle F 2 . The bisecting line £ between 
0 and V in the plane of F 2 contains p. The line £ intersects two edges of F 2 , once 
at v/2 and once at some other point p'. If p is (strictly) outside F 2 , then |p| > \p'\. 
This leads to a contradiction, once we show \p'\ > y^. If p' lies on the edge between 
V and W 2 , this is clear, because an elementary exercise shows that every point on 
the line passing through W 2 and v has distance at least from the origin. Assume 
p' lies between 0 and W 2 , and consider p' as a function of W 2 and v. Its length \p'\ 
attains its minimum when F 2 is the right triangle |n| = 2 , \v — W 2 \ = 2 , 1 ^ 2 ! = 2 v^. 
Thus, Ip'I > v^. 
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So p lies in the interior of F 2 . The vertex wi has distance at most \/2 from p. No 
vertex wi can come within \/2 of an interior point of F 2 nnless the circnmradins of 
F 2 is at least v^. If the circnmradins of F 2 is then p is the circnmcenter of S', so 
that the circnmradins of S is v/2, making S a small simplex. If the circnmradins is 
greater than then since \p — wi \ < \/2, p lies in the Voronoi cell at wi. Thns, wi 
clips (possibly degenerately) a standard region across the faces F 2 from rci, based 
at 0, n, or W 2 - By Lemma 2.2 (rci, 0, n, rc 2 ) is a qnasi-regnlar tetrahedron. I 


We continne with onr description of the hgnre in Diagram 4.3. The arc a cannot 
be interrnpted by a farther (Voronoi) edge of F'. Snch an edge wonld be dnal to a 
(Delannay) face with vertices 0, n, and some v'. The circnmradins of the triangle 
with these three vertices wonld be less than (becanse every edge in F' comes 
within distance y^ of the origin). This contradicts the constrnction of F" with 
half-planes given above. This completes onr discnssion of the hgnre in Diagram 
4.3. We emphasize, however, that the edges ei or 62 may degenerate to length 0, 
and the circnlar arc a may degenerate to a point. 

This Voronoi face-wedge can be broken into three convex pieces: the convex hnll 
of 0, n/2, and the circular arc, and the convex hulls of 0, n/2, and the edge e^, 
for i = 1,2. The hrst piece has the third type, the others have the fourth type. 
The boundary condition 4/3 < expresses the fact that the circnmradins of a 
triangle with sides of length at least 2 cannot be less than 2y^/3. This completes 
the reduction to the four given types. 

Now we must show that each of the given types has density at most 5oct- 

Type 1 : A small simplex that is not a quasi-regular tetrahedron. Let S' be a 
small simplex of Voronoi type with at least one edge longer than 2.51. By the 
monotonicity properties of the circnmradins, we know that the circnmradins of S 
is at least rad(5'(2, 2, 2, 2.51, 2, 2)) > 1.3045. Let 5(a, 6 , c) denote the density of the 
Rogers simplex R(a, 6 , c) (see 1.8.6). By Rogers’s lemma (1.8. 6 .2), the six Rogers 
simplices Vs have density less than 5oct and vor(5', V 5 ) < 0 if the circnmradins of 
the three faces is at least 1.207 (h(l, 1.207,1.3045) < 5oct)- This condition on the 
circnmradins of the faces holds whenever there are two edges longer than 2.51 at 
the origin ( 77 ( 2 .51, 2, 2) > 1.207) or whenever there are two oppositely arranged 
edges longer than 2.51. 

Thus, to show that vor(5') < 0 for small simplices of Voronoi type, we must consider 
the following cases: (1) one edge longer than 2.51, (2) two adjacent edges longer 
than 2.51, and (3) three edges longer than 2.51 meeting at a vertex. These cases 
are covered by Calculation 4.5.2. In (1), we may assume that at least one of 
the conditions for compression scoring in Section 3.1 fails to hold. In Calculation 
4.5.2.2, we may make the stronger assumption rad(S') < 1.39, for otherwise, the 
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Rogers simplices at the origin have density at most h(l, 17 ( 2 , 2, 2.06), 1.39) < doct so 
that vor(5') < 0 . 

We rely on Calcnlation 4.5.1 for small simplices of compression type. The appendix 
proves the result for simplices in an explicit neighborhood of 5'(2v^, 2 , 2 , 2 , 2 , 2 ). 
These calculations are established by methods of interval arithmetic described in 
Part I. Source code appears in [ 6 ]. 

Calculation 4.5.1. If S' is a small simplex that is not a quasi-regular tetrahe¬ 
dron, then r(S) < 0. If equality is attained, then the simplex S is congruent to 
S( 2 v^, 2 , 2 , 2 , 2 , 2 ) or to the simplex of zero volume S( 2 v^, 2 , 2 , 2 a/ 2 , 2 , 2 ). 

Calculation 4.5.2. Assume S is small. vor(S(yi,... , 7 / 6 )) < 0 if (7/1,... , 7 / 6 ) 
belongs to any of the cells (l)-(ll). Let I denote the interval [2,2.51] and L = 
[2.51,2^2]. 

(1) L[2.06,2.51]/4, 

(2) L/2[2.06,2.51]/2, 

(3) L/3[2.08,2.51]/, 

(4) [2.06,2.51]/2L/2, 

(5) /[2.08,2.51]/L/2, 

(6) /3L[2.2,2.51]/, 

(7) /3[2.58, 2^2] [2.12, 2.51]/, 

( 8 ) LPLI, 

(9) 

( 10 ) PL‘^1, 

(11) /^lA 

Type 2 : The set {tx :0<t<l, x E P 2 }- In this case the density is increased by 
intersecting the set with a ball of radius 1.18 centered at the origin. The resulting 
intersection has density 1/1.18^ < doct, as required. 

Type 3 : The set {tx :0<t<l, x E P 3 }. The bounding circular arc of P 3 has 
distance from the origin. The set has the same density as a right circular cone, 
with base a disk of radius \/2 — and height h. This cone has volume 7 i{2 — h‘^)h/3. 
The solid angle at the apex of the cone is 27r(l — cos 6*), where cos 6* = hj^/2. This 
gives a density of \/2/{h‘^ -|- h\/2). This function is maximized over the interval 
[1,1.18] at h = 1. The density is then at most 2 — ^J2 < 5oct- 

Type 4: A Rogers simplex Ri = R{a, 6, y/2). where 1 < a < 1.18 and 4/3 < 6^ < 2. 

By Lemma 1.8.6.2, the density of this simplex is at most that of the Rogers simplex 
/?2 = /?(l,2v^/3, V 2 ). This simplex has the density Soct of a regular octahedron. 
(In fact, the regular octahedron may be partitioned into simplices congruent to R 2 
and its mirror.) We see that the original simplex Ri has density 6 oct exactly when. 
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in the notation of 1.8. 6 .2, |si| = |s 2 |, for all Ai, A 2 , and A 3 as above. This implies 
that a = 1 and h = 2^/3/3. This completes the proof of Theorem 4.1. I 

Proposition 4.6. A clnster other than a qnasi-regular tetrahedron attains a score 
of Qpt if and only if it is made up of simplices congruent to 5'(2,2,2,2,2, 2 \/ 2 ), and 
possibly some additional simplices of zero volume. 

Proof. Types 2 and 3 always give strictly negative scores for regions of positive 
volume. According to Calculation 4.5, a region of the hrst type with positive volume 
gives a strictly negative score unless it is congruent to S{2, 2, 2, 2, 2, 2\/2). 

Consider a region of the fourth type with score Opt. We must have a = 1 and 
b = 2a/3/3. The circumradius of the faces Fi and F 2 of Diagram 4.3 is then 2-\/3/3. 
This forces the faces Fi and F 2 to be equilateral triangles of edge length 2. The 
arc a in Diagram 4.3 must reduce to a point. The edges ei and 62 in Diagram 4.3 
“ if they have positive length - must then meet at a point at distance ^/2 from the 
origin. This point is a vertex of a Voronoi cell and the circumcenter of a Delaunay 
simplex S (of circumradius a/2). The only simplex with two equilateral faces of side 
2 and rad(5') = y/2 is the wedge of an octahedron S = S{2, 2, 2, 2, 2, 2a/2). This is 
a small simplex. 

The other possibility is that both the arc a and an edge (say 62 ) degenerate to 
length 0. In this case, Lemma 4.4 shows that the restricted cell belongs to a small 
Delaunay simplex or a quasi-regular tetrahedron. These cases have already been 
treated. I 
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Appendix 


We give a direct argument that r(5') < Opt, when the lengths of a small simplex S 
are within 0.001 of S'o = S{2^/2, 2, 2, 2, 2, 2). Set S = S{yi, 1 / 2 , ys, 1 / 4 , 1 / 5 , l/e)- Write 
yi = 2\f2 — f\, and yi = 2 + fi, for i > 1, where 0 < fi < 0.001. Set xi = y^ = 8 —ei 
and Xi — yf = 4 + Ci, for f > 1. Then 0 < < 0.006. Recall from Section 1.8.4 

that 


1 2 , „.2 „. 2 \ I t „.2 I „.2 „. 2 \ I ,'„.2 I „.2 


a( 2 /i, 1/2, ■ • ■ , l/e) = 1/11/21/3 + 771/1(1/2 +1/3 -1/4) + 771/2(1/1 +1/3 -1/5) + 771/3(1/1 +1/2 -l/e)- 


Set 


ao = a(|/i, 1/2, 1 / 3 , 1 / 4 , 1 / 5 , l/e), 
= a(l/l, 1 / 5 , 1 / 6 , 1 / 4 , 1 / 2 ,1/3), 
02 = a(l/ 4 , 1 / 2 , 1 / 6 , 1 / 1 , 1 / 5 , 1/3), 
O3 = 0(1/4, 1 / 5 , 1 / 3 , 1 / 1 , 1 / 2 , l/e). 


aoo = 0(272,2,2,2,2,2) 

oio = 0(272, 2, 2, 2, 2, 2) 
a 2 o = 0 ( 2 , 2, 2, 2^2,2) 

aao = o(2,2,2,2v^,2,2) 


16 + 1272, 
16 + 1272, 
16, 

16. 


Section 1.8.4 and the bounds on fi give ai > a- , where = 32.27 and 

a 2 = O 3 = 15.3. 

Let A be the function of Section 1.8.1, and set Aq = A( 8 , 4,4, 4,4, 4). Set t = 
7A(a:i, ... , xq)/2, and to = 7^0/2. A simple calculus exercise shows that 

A(a:i,... , xq) > A( 8 , 4,4, 2 : 4 , 4,4) = 128 — 8 e|. 


This gives t > 5.628. Let bi = 2/(3(1 + to/^io))’ 60 = 61 = (3 + 272)/9 and 

b 2 = bs = 16/27. Set cq = —Soct /0 + ~ —0.00679271. Then 


(t - to)co 


(A - Ao)co ^ -2e|co 
4(t + to) {t + to) 


< 0.002e^ < 0 . 0002 / 4 . 


We are ready to estimate r(S'). An argument parallel to that of Lemma 1.9. 1.1 
gives 


(1) ns) < ns„) + (t- to)co+<7 °‘* + * 1 : ~ 

The two sums on the right-hand side are polynomials in fi with no constant terms. 
To give an upper bound on these polynomials, write them as a sum of monomials, 
and discard the negative monomials of order greater than 2. The positive monomials 
of order greater than 2 are dominated by 

/7/2^ ■ ■ ■ ft < {omi)<^^+-+<^^-\h + /2 + ■ ■ ■ + /e). 

This approximation shows that the hrst sum in (1) is at most —0.005/i — O.O 4/4 — 
0-03(/2 + /3 + 7 + /6) and the second sum in (1) is at most 0.00056(/i + /2 + - ■ ■T/e)- 
The result easily follows. 

This argument is easily adapted to a neighborhood of Si = 5'(272, 2, 2, 272, 2, 2). 
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Co ~ —0.0225, and t > 0. A similar argument leads to the conclusion that r(5') < 
r(5'i) = Opt, if A is a small simplex such that S Si, and the lengths of the edges 
of S are within 0.01 of those of Si. 
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